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ABSTRACT 

Skewness of temperature fluctuations of the cosmic microwave 
background (CMB) produced by initially Gaussian adiabatic 
perturbations with the flat (Harrison-Zeldovich) spectrum, which 
arises due to non-linear corrections to a gravitational potential at the 
matter-dominated stage, is calculated quantitatively. For the standard 
CDM model, the effect appears to be smaller than expected previously 
and lies below the cosmic variance limit even for small angles. The sign of 
the skewness is opposite to that of the skewness of density perturbations. 



Subject headings: cosmology: theory — cosmic microwave background 
— large-scale structure of the Universe 
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1. Introduction 



The reliable detection of fluctuations of the CMB temperature AT(9,ip)/T 
at large angular scales by the COBE group (with the data being in a very good 
agreement with a prediction made on the basis of the inflationary scenario of the early 
Universe 10 years before the detection) stimulates further investigation of subtler 
effects. Among the most important of them are possible deviations of the statistics 
of these fluctuations from the Gaussian one. The basic result which follows from 
all sufficiently simple variants of the inflationary scenario (contrary to that of rival 
theories based on topological defects such as cosmic strings, etc.) is that the statistics 
of the AT/T fluctuations is Gaussian because they are linearly connected to quantum 
vacuum fluctuations of a very weakly interacting scalar field (the inflaton). Thus, in 
the leading (linear) approximation, the mean CMB skewness 6*3(0) = ((AT/T) 3 ) = 
where () denotes averaging with respect to different realizations of stochastic 
space-time metric perturbations of the Friedmann-Robertson- Walker cosmological 
model which produce AT/T. Note that the corresponding observable quantity is 
C 3 (0) = (An)' 1 f(AT(6,tp)/T) 3 dQ, i.e., the average over the sky of one particular 
realization. Therefore, C 3 (Q) generally differs from C 3 (0) by a so-called cosmic 
variance which is of the order of (((AT/T) 2 )) 3/2 (see, e.g., Srednicki 1993). In 
practice, even 6*3(0) is unachievable because of a finite beam width of antennas, 
incomplete sky coverage, etc. 

There exist different physical effects which may result in the appearance of a 
small, but non-zero mean skewness. Some of them are connected with non-linear 
corrections to the initial spectrum of scalar (adiabatic) metric perturbations 
which were generated during inflation (Falk et al. 1993, Gangui et al. 1994), the 
corresponding part of the CMB mean skewness may be called primordial. Other 
effects which take place after recombination (Luo & Schramm 1993) produce a 
secondary skewness. In this paper, we will concentrate on a detailed calculation 
of the secondary skewness produced by non-linear corrections to the primordial 
gravitational potential <3> which arises due to the same gravitational instability at 
the matter-dominated stage that leads to formation of galaxies and the large-scale 
structure of the Universe. The corresponding contribution to AT/T is contained in 
the non-local term of the Sachs- Wolfe expression for AT/T (Sachs & Wolfe 1967), 
it is also called the Rees-Sciama effect (Rees & Sciama 1968, a different view on this 
effect is presented in Zeldovich & Sazhin 1987). 

The reason for our primary interest in this effect is that its contribution to 
AT/T, though formally being of second order in powers of a small initial gravitational 
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potential $(0) = <po{r), is not much less than the main linear effect. Its value may 
be estimated as AT/T ~ <Mp/p where Sp/p is the present rms density perturbation 
inside some linear characteristic size L. For L ~ R eq rs 30/i _1 Mpc corresponding to 
the angular scale $ ~ LH /c ~ 30' (if viewed from the redshift z ~ 1), it is ~ 0o-2e g 
where z e(? ~ 4- lO 4 ^ 2 /?" 1 is the redshift of matter-radiation energy equality; therefore, 
it may reach 10~ 6 (Martinez-Gonzales, Sanz & Silk 1992). Here h is the value of the 
Hubble constant H in units of 100 km/s/Mpc and k = e m d/e 7 = 1.68 for 3 types of 
light (m < 1 eV) neutrinos with standard concentrations. 

On the basis of the above argument, it was conjectured (Luo & Schramm 1993) 
that the secondary CMB skewness produced by the non-local part of the Sachs- Wolfe 
effect will dominate the primordial skewness generated in the inflationary scenario. 
However, our calculations presented in the next section show that this is not the case 
for the standard CDM model. The secondary skewness appears to be small and of 
the order of primordial skewness. A possibility of getting a larger skewness without 
introducing late-time phase transitions is discussed in the Sec. 3. The section also 
contains our conclusions, as well as comparison with the cosmic variance of the 
skewness in the case of purely Gaussian perturbations. Details of our calculations 
are displayed in Appendices A and B. 



2. Method of calculation 

For large angular scales $ > 2°, CMB temperature fluctuations produced by 
adiabatic perturbations are given by the Sachs- Wolfe formula: 

V T J loc \ T ) non-loc 3 J Vrec \ Of] ) r=m ^ 

(r,9,ip), r = 7] - 7] rec , r/ = y"-^y, (1) 

where a(t) is the scale factor of the FRW model, r] is the present conformal time 
and r] rec is the conformal recombination time, T] /T] rec w z)H ~ 30. More exact 
calculation for the CDM model with 3 light neutrino species and z rec = 1100 
gives f]o/Vrec = 49.6. At the matter-dominated stage in the absence of spatial 
curvature and the cosmological constant, a(t) oc t 2//3 , 77 = 3t/a(t). For smaller angles 
5' < •& < 2° where standard recombination may be still considered as instantaneous, 
additional local terms in Eq. ([]]) appear which describe the Doppler and the Silk 
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effects (see, e.g., Starobinsky 1988), but the non-local term remains the same. In 
the linear approximation, this term is exactly zero for any power-law a(t), though 
it may contribute significantly if a(t) deviates from a power-law behaviour after 
recombination, e.g., due to the cosmological constant (Kofman & Starobinsky 
1985), decaying relativistic particles (Kofman, Pogosyan & Starobinsky 1986) 
or spatial curvature (Wilson 1983, Abbott & Schaeffer 1986). However, for the 
purely matter-dominated stage in the flat Universe that we consider, the non-local 
term produces a non-zero contribution to AT/T if non-linear corrections to the 
gravitational potential <£>, arising due to gravitational instability in the Universe (see, 
e.g., Peebles 1980), are taken into account. These corrections are also responsible 
for the appearance of a non-zero skewness of CMB in the case of a Gaussian initial 
(linear) potential. 

So, we expand the (peculiar) gravitational potential into a series in powers of 
a density enhancement: $(r, t) = + & 2 ' + .... The linear term $W = 0o( r ) is 
assumed to be Gaussian. Further, we consider the case of a flat (Harrison-Zeldovich) 
initial spectrum. Then the two-point correlation function of the linear potential is 
given by 

Mr) = (M0)Mr)) = fp,(kf^k*dk = Bj c 2 (k)^f , (2) 

where c(k) is the standard transfer function of the CDM model. For physical scales 
R = ra(t ) ^> R eq , £^(r) ~ £^(0) — B\n(r/r eq ) (strictly speaking, the constant 
£^(0) is infinite but it is unobservable, so no difficulties arise). The constant B is 
related to the amplitude normalization A introduced in Starobinsky 1983 by the 
relation B = 9 A 2 /2007T 2 . The rms CMB quadrupole value is expressed through it by 
Qrms-Ps/T 2 = 55/108, so that B = 1.16 x 10~ 9 for the presently preferred values 
T = 2.726 K and Q rms ~ps = 20 \iK (see, e.g., Gorski et al. 1994). 

The second-order term is given by (Peebles 1980, notations of the paper by 
Munshi & Starobinsky 1994 are used below): 

P(r) = (A0 O ) 2 + V0oV(A0 o ) = V(V0 o A0 o ), Qif) = (V0 O ) 2 - (3) 
To get a finite expression for C 3 (0), it is sufficient to use the "renormalized" quantity: 

^ n (r,v) = ® {2) (r,v)-(® {2) (r,v)) , (4) 

(the last term in (PS) depends on 77 only) that corresponds to the subtraction of an 
unobservable constant (monopole) term from AT/T. More rigorous justification of 
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this prescription follows from the fact that the difference $(r, 77) — $(r ,^), where r 
denotes an observer (e.g., our) location, is an observable and finite quantity both 
in linear and non-linear regimes. Further, we omit the subscript "ren". Note that 
the subtraction @ does not remove the whole contribution to 6*3(0) from monopole 
terms, in each of AT /T some finite part remains. 

A non-zero contribution to the CMB mean skewness in the lowest (fourth) order 
in $ is given by 



Let us now remove a monopole component from each of 4£ completely: 



loc \ J- / non—loc 

\ 

T 

AT\ AT 1 r AT 
~Y)s ~ J ~T~ 



dQ . (6) 



Then the term (|5]) can be represented as 

^3^(0) = luuu — Iuum — %Iumu + 2Iumm + Immu — Immm (7) 

where the subscript U means substitution of an unsubtracted into (|) and M - 
substitution of a monopole term there (the gravitational potential in both terms is 
assumed to be renormalized according to (|J)). Strictly speaking, a dipole component 
should be subtracted from further, but this correction appears to be smaller 

than the correction due to the monopole subtraction @ (e.g., Iuud ~ 0.74 Iuum, 
see Appendix A) and practically does not change the final result for 6*3(0). 

The main term in the right-hand side of Eq. ([/]) is luuu- After substitution of 
(H) into Eq. (|l|), it takes the form 

W = ^/ ^fe 2 + 5A- 1 ((A^) 2 + ^(A^)') dr, (8) 

where ^ is given in Eq. (0) and the prime means derivative with respect to 
r. Note that for r > r eq , A" 1 ((A^) 2 + ^(A^)') = -Jj. Thus, if we are 
speaking about the skewness of temperature fluctuations smoothed over a scale 
R s > R rec = a(r] )r] rec (e.g., due to a finite antenna beam width), then the integral in 
(|8|) diverges at both large and small r\ logarithmically: luuu — ~ ^"-^ 2m 'i k ) where 
Rh = a(r] )r]Q is a present cosmological horizon of the FRW model. This means that 
the growth of ( 4£) in (B) at late times due to <3>( 2 ) is exactly cancelled by 

V 1 / non—loc 

decay of correlations between (4^ ) and (^-) at large spatial separations. 

V 1 / loc \ 1 J non—loc 

Another conclusion following from the behaviour of the integrand in Eq.(^) is that 
the contribution of the small redshift region z ~ 1, rj ~ rjo to the skewness is 
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sub dominant, it is smaller than the total effect in a few (~ 1/ \n(r] /r] rec )) times. 
That is why higher-order non-linear corrections do not noticeably change the fourth 
order result for the skewness (we check it explicitly for the leading sixth order term, 
see Appendix B). Clearly, the same statement is true for the case of a blue-tilted 
initial spectrum n > 1. However, in the opposite case of a significantly red-tilted 
spectrum (n ~ 0.9) higher-order non-linear corrections may become crucial since the 
main effect originates at recent times. 



In terms of the multipole decomposition, this behaviour corresponds to equal 
contributions to 6*3(0) from each logarithmic interval of I for I ^> 1. In other words, 
if 6*3(0) = Y,iCi, then Ci = — 2 i(i+o 5) f° r ^ ^ 1- Using the monopole value 60 



obtained below (see Appendix A and Eq. fllTD), this fit may be made better by 
changing (Z + 0.5) to (/ + 1/3). So, if we take the smoothing angle FWHP = 10° as 
in the smoothed COBE maps [6 g = 0A259fwhp — 1/13.5), then 

63(10°) = -0.165 2 , ((AT(10°)) 3 ) = 4.4 (^K) 3 (9) 



(with monopole and dipole terms subtracted). Since 6 2 (10°) = 0.185, the smoothed 
large- angle skewness parameter S3 (10°) = ^|^ j ~ — 5. 

For the unsmoothed distribution of temperature fluctuations, the main 
contribution to Ijjuu comes from the vicinity of the recombination surface 
V — Vrec ~ Veq as a result of properties of the transfer function c(k). However, this 
distance is still significantly larger than the thickness of the recombination surface. 
Hence, the latter may be neglected for our problem. The corresponding angular 
range is 10' — 30'. Here lies the first acoustic, or Doppler peak, so the first term alone 
in the right-hand side of Eq. ([!]) should be corrected by accounting for the Doppler 
and the Silk effects (but not the second one). As a result, we obtain (details of the 
calculation are presented in the Appendix A): 

Iuuu ~ -2.2 B 2 . (10) 

Other terms in the right-hand side of (0) may be calculated in the approximation 
r] rec = 0, c(k) = 1. From symmetry considerations, Ijjmm = Immu = hiMM- It 
is shown in the Appendix A that I UUM = -^(ir _ ln2)£? 2 -0.055 2 . Other 
auxiliary integrals entering into Eq. (|7|) are also given there. The dipole contribution 
can be estimated to be ~ — 0.07B 2 using the abovementioned fit Ci = — 21(1+1/3) ■ ^ s 
a result, 



63^(0) _ Iuuu ~ Iuum — 2(Iumu ~ Iumm) — 61 



-7- 



= -2.2B 2 + 0.055 2 + 0.245 2 + 0.07E 2 « -1.8B 2 , 
((AT) 3 ) « 50(/^) 3 . (11) 

With a good accuracy, <j\ = C 2 (0) = ((^) 2 > « lO" 9 « B where 
C 2 ('#) = (^r-(O) ^-(i?)) is the 2-point angular temperature correlation function. 
Therefore, the skewness parameter S* 3 (0) = ~ — 2 is of the same order as the 
primordial skewness considered in Gangui et al. 1994. The unsmoothed value of the 
skewness 6*3(0) is ~ 11 times more than 63(10°) because the main effect comes from 
the angles 10' — 30'. However, the unsmoothed skewness parameter $3(0) is smaller 
than 6*3(10°) because fluctuations themselves (i.e. C%) are significantly larger at 
small angles. 

Note the negative sign of S3. Its physical explanation is that it reflects the 
existence of large regions with $ positive and growing with time (corresponding to 
voids) which produce positive (4?-) and relatively small regions of negative 

V ' / non—loc 

and decreasing $ which produce smaller cold spots but with larger absolute values 
of 4rr. It is these cold spots that make the main contribution to the skewness 
and determine its sign. The sign is opposite to the sign of a skewness of density 
perturbations ^. In the latter case, the main contribution to the skewness is 
produced by regions with & > which are smaller in volume but larger in amplitude 
of 1^1. 



As noted above, the smallness of |5*3| is due to the fact that the region where 
l&^l along a photon trajectory becomes comparable to |0o| [z ~ 1, f] ~ ^0) is 
widely spatially separated from the recombination surface 11 = ri rec where (4?) is 

V 1 J loc 

located. Thus, to check that there is no other significant contribution to |5*3|, we 
consider sixth-order terms that need not possess this property. The most important 
of them is the following term: 

^ 6) (0) = ((^) 3 ), (12) 

\ -l / non—loc 

where (4?) is calculated with the use of $^ 2 \ The main contribution to (1X21) 

V T Jnon-loc V| v 

is mainly produced at recent times 77 ~ rjo and it is not attenuated by a small value 
of the spatial correlation function. However, it contains one more power of B as 
compared to (|TTD. The term ([12]) can be represented as 

g3 6) (0) = - 1 fnWn 8 N7 / / / ^ hd* k 2 d 3 k 3 (h )<f (h)^ ( k 3 )6(k lz - k 2z )6(k 2z ~ k 3z ) 



42 3 (2tt)' 



xM(k 1; -k 2 )M(k 2 , -k 3 )M(k 3 , -k!)(13) 
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where the kernel M(k 1; k 2 ) is defined in Eq. (|A8|) below (see Appendix B for a 
detailed derivation). Numerical evalulation of this integral gives 

Cf )(0) = -2.15 x 10 6 B\ Sf 1 ps -lO- 3 . (14) 

Thus, the sixth order contribution is much smaller than the fourth order one. 



3. Conclusions 

We have calculated the mean CMB skewness generated due to leading non-linear 
corrections to an initially Gaussian adiabatic perturbations with the flat spectrum. 
The unsmoothed value S^O) ~ —2 that we found lies much below the cosmic 
variance of this quantity 5 S3 ps (cWc) 1 ~ 130 (here we model C^t?) by the Gaussian 
C 2 = <rf.exp (—^p-) with ot = 3 • 10~ 5 and l c = 250 corresponding to the first 
Doppler peak). The situation for smoothed large-angle maps is even worse (e.g., 
S 3 ps —5 but 5 S3 ps 3600 for 9fwhp = 10°). Therefore, as regarding observations, 
the prediction is that there should be no noticeable mean skewness above the noise 
level due to cosmic variance in the case of the standard CDM model with the flat 
(Harrison-Zeldovich) initial spectrum of adiabatic perturbations. This conclusion 
seems to be in a good agreement with existing data ( Hinshaw et al. 1994, Kogut 
et al. 1994). It may be considered as one more confirmation of predictions of the 
inflationary scenario, though, of course, this fact does not close the way for other 
theories leading to the same prediction. 

It is clear from our derivation that IS3I can be substantially larger if there 
exists a first order contribution to at late times (z ~ 1), because then 

\ ' / non—loc 

correlations between first and second order terms in 6*3(0) are not small. It may 
happen, as mentioned above, for the flat CDM+A cosmological model and other 
ones. We shall consider this case elsewhere. 
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A. 



Fourth order contribution 



We outline some of the details of our estimation of the skewness. The lowest 
order term leading to skewness (4th. order in </> ) is given by 



where 0o an d ( 3^ 2 ' ) represent values of potential fluctuations at the linear and 
second order, respectively, and <3>( 2 ) is regularized according to Eq. (£|); r is 
defined in Eq. ([!]). In the paper, we follow the notation 4> = 4>u ~ 4>m — 4>d and 
$( 2 ) = <s> v — $ M — where the subscripts U, M and D represent unsubtracted, 
monopole and dipole fields, correspondingly. 

The CMB mean skewness at the lowest order can be then expressed (see Eq. 



where we have ignored terms involving <pr> since they are expected to be small. In 
addition, we also expect that Iumd and Immd are smaller in magnitude. 

All the terms contain an average of the type (0o( r i)0o( r 'i) ( ^ )( ' 2 ^( r 2)) which can 
be expressed in terms of the potential-potential correlation function £^ and its first 
derivative £i as 



(0 o (r 1 )0o(r' 1 )$ (2) (r 2 )) = 2V^(r x - r 2 )V^K - r 2 ) + [ 2^(n - r 3 )^(r[ - r 3 ) 

+V£»(r 1 - r3)Ve;(r; - r 3 ) + V&fa - r 8 )Vfc(ri - r 3 ) ] . 



The terms Iuum and Iuud involve a simpler expression with ri = r^ and can be 
easily evaluated in the coordinate space. We find that the terms Iuum, Iuud, Iumu 
and Iumm are all insensitive to the form of the transfer function c(k) and may be 
evaluated in the approximation c(k) = 1, i] rec = 0. 




(Al) 



(0)) ^ 
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The term Ijjum and Ijjud involve the monpole and dipole term of &( 2 \ 
respectively and can be expressed as 

2 1 /■ rna , d 



lUUM — 



63 4tt 



't/t/D 



dn 2 / ^^(^( ri )$ (2) (r 2 )) 

Jrjrec (J 1 1 



Assuming a scale invariant spectrum and c(/c) = 1, 77, 
£^(r) ~ £^(0) — -Bln(r/r e9 ) to reduce the expressions to 

2B 2 1 /■»» , r dVL 

IUUM = —— J o d V V J ~2 , 

cos 6 



|ri|=r , \r 2 \=ilo-V ' 
|ri|=n>, |r 2 |=r?o-r? 

= 0, we set 



(A3) 



IjJUD — 



21 4tt 
2B 2 3 /■ cos 

■y ^w- 



21 4tt 



(A4) 



where r 2 = (770 — rj) 2 + rj 2 — 2r] (r] — rj) cos 6*. 

Performing the integration over angles, we get 



WUM 



B r 1 x dx , (2 — x 
In 

V x 



21 Jo 1 - x 



25 2 - 2 



,7T 



I ~ 52 



1 xdx 



7 Jo (1-x) 



'l-x) 2 + l , (2-x 
In 



x 



1 +x 



( In 2) « -0.05155 2 , 

21 v 8 ; 

"T [I " t) w -°- 0380i?2 ( A5 ) 



Now we consider the terms ij/Mt/ and Iumm which involve the monopole term 
of 0o and can be written as 

2 1 r P° , 9 



WMU 



63 4tt 



^ / dry— (MnMriWra)) > 

Vo-V 



|n| = ir^ = r , r 2 = ri 

»7o , d 



r 



L UMM 



2 \ c c r 7 !® n 

63 IS? I dS1 '' l' m \L *%(Mr,)Mr[)^(r 2) ) , 



n = r x = r , r 2 = Vo ~ V ■ 



(A6) 



It is more convenient to evaluate these terms in the momentum space. In the k - 
space, the unsubtracted, monopole and dipole terms for the temperature fluctuation 
along the line of sight n can be expressed as 



m: - y 



d^k 

J2nf 



We 



jknro 
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AT V 1} 1 ! * h * (v\ ■ (U 1 
-T)m = 3 7 (2^ 0o(k)jo(A;ro) 



/AT\( 2 ) 2 r d 3 k , 9W , i rvo 



where 



\ w / r (1 k rvo 

Vr')„ = 21 /(2^* <2,(k) X.„ *>■>*<*<»->)>■ 

/AT\( 2) 2z /• rf 3 A; ,, 2Wl . kn f?o . , .. 

{Tin = 7i^ ()(k) TL^^-^ (A7) 

$(2) (k) = - ?*\\?* k * M(k 1; k 2 ) ^(kOMka), 



(2tt) 3 v 7 (27r) 3 (27r) 3 

M(ki,k 2 ) = 5P(ki,k 2 ) + Q(ki,k 2 

2k\k 2 + k^k 2 + k? 
2(k! + k 2 )2 



^(ki, k 2 ) = 1 2 J, V 2 ; * ~ » Q( k ^ k 2) = k ik 2 , k = \k\, (AS 



the function j n (x) represents the standard spherical Bessel function of order 
n. Note that within the scope of this paper the first order temperature 
fluctuation (AT/T)^ = (AT/T)/ oc and the second order temperature fluctuation 
(AT/Tp 2 ^ = (AT/T) nori _ loc , owing to the fact that the linear gravitational potential 
does not change in a dust-dominated, flat FRW universe. 

In evaluating expectation values in the fc-space we invoke the Gaussian nature 
of the initial potential perturbations 4>o and implement the following relation 



<>o(kO0o(k 2 )<I> (2) (k)} - (0 o (k 1 )0o(k 2 )0o(k 3 )0o(k4)) = (j) 2 (ki)4> 2 (k 
x (5(k! + k 3 ) 5{k 2 + k 4 ) + 5(k x + k 4 ) 5{k 2 + k 3 )) , 



2ir 2 Bc 2 (k) 

<P\k) = k3 (A9) 

where we have ignored the self-coupling term <p 2 (ki) 5{k\ + k 2 )($( 2 )(k)) because it 
is just cancelled after the renormalization (^). 

The expressions for the two terms now read 



l UMU 



b 2 r dr]V f d ^i fd^^Mi ^) M(kl)k2)e -^ in (,-^) elk2n( , -,) io( ^ o) 



63 Jrjrec J 2tc J 2tt k\ k\ 
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We find it convenient to split each of the terms into two pieces involving the 
local part Q and the non-local part P part of $^ 2 - ) , respectively. After carrying out 
some of the integrations analytically, the terms can be expressed as 



AB 2 rvo 



*-UMU Q 



63 



rrjo roc roc 

/ di]ri / dki c (h) ji{ki{r]-r] rec )) / dk 2 c (k 2 ) jo(k 2 r ) ji(k 2 {r]o-rj)) , 

J ri rec JO JO 



5B 2 no 



WMUp — 



63 



x J' du \( k l + k2 2) u + 2k ^ 



k\ + k 2 + 2k 1 k 2 u 



I drj r] / dki c 2 (ki) \ dk 2 c 2 (k 2 ) j (k 2 r ) 

J r) rec JO JO 

jo [\lkl{r]- r) rec ) 2 + kl(r]Q - rf) 2 - 2k 1 k 2 {r] - rj rec )(r]o - rf) 



u 



4B 2 fvo 

63 Jri rec 



Iummq = I drj r] 



dk c 2 {k) j (kr ) ji(k(r} - r})) 



UMMp 
1 

X / (III 



5B 2 f^ f f 

= — ^- / dr] r] dkic 2 (k l )jo(k l ro) \ dk 2 c 2 (k 2 )j (k 2 r ) 

UO Jr) rec J J 



The expressions are further simplified if we take c{k) = 1 and r\ rec = to obtain 
the final results 



2B 2 r 1 dx 



l UMU Q 



63 



r- 

Jo 1 



X 



„ x(x — 2) , ! . 
1+-^ tanh _1 (l-x) 



2B 2 



63 



(21n2-l) « -0.0123B 2 , 



Iumup — — 0.23-B 2 , J, 



UMU — 1 UMU Q + IuMUp 



-0.2AB 2 



(A12) 



and 



5 5 



'C/MMq 



1 £ 



63 Jo (1 - xf 



1 + 



— 2) 



tanh *(1 — x) 



7.30 x 1(T 4 £ 2 , 



luMMp ~ — 0.125 2 , Iumm — Iumm q + Iumm p — — 0.12B 2 . (A13) 
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We have numerically verified the above terms including the CDM transfer function 
and value of r\ rec . 

Now we deal with the most significant term Iuuu which depends sensitively 
on small-scale power in the radiation perturbation spectrum and transfer functions 
for radiation (see below). As noted before, for the angular range 5' < d < 2° that 
roughly corresponds to the scale range 10 Mpc < a^^i^kh)^ 1 < 200 Mpc, all three 
effects - the Sachs- Wolfe, the Doppler and the Silk ones - should be taken into 
account in the local term (AT/T) /oc . Moreover, only the sum of the three effects 
appears to be a gauge-invariant quantity. On the other hand, the recombination may 
be still considered as instanteneous (the width of the recombination surface in terms 
of the conformal time Ai] rec <C A; -1 ). More exactly, the effect of the non-zero width 
Ai] rec may be empirically accounted by the increase of the scale of Silk damping. 
Therefore, we use (partly unpublished) results of Starobinsky & Sahni 1984 and 
Sahni 1984 (see also Starobinsky 1987, 1988) obtained in the two-fluid approximation 
of the CDM model with radiation. Namely, the limiting case £7^ = was first 
considered and the matter before recombination was assumed to consist of two 
ideal fluids interacting through gravity only: dust with pressure p = representing 
cold dark matter and radiation with p = e/3 representing photons (tightly coupled 
with baryons) and other massless particles. After the recombination, photons are 
described as free massless particles. This approach is equivalent to that used in 
Seljak 1994. Also, we assume that r] eq <C r] rec . However, a small, but non-zero value 
of Qb should be finally taken into account because it is crucial for the determination 
of the Silk damping scale. Also, we accounted for it in the value of the sound velocity 
in the radiation component before recombination to get the right location of the 
acoustic peaks. 

We use the following values in actual caculations: Hq = 50 km/s/Mpc, 
T = 2.726 K, k = 1.681, Q tot = 1, Q b = 0.06, z rec = 1100. The scale factor of the 
two-fluid (dust and radiation) CDM model has the form aijf) = 0177(77 + rji) where 
r]i = 2(y / 2 + l)r] eq = A.828r/ eq . Now we can find values of r/ eq , r/i and r\ rec using the 
present-day value of Q 7 = 0.9895 x 10" 4 : 

— « (2^K i Y l (l - y^/) = 38.3 , — = 185 , z eq = (kQ 7 ) _1 = 6012 , 

— = ¥ — — = 0.771 , — = 49.6 . A14 

m 2 r] rec 

It is clear that the approximation r] eq 77 rec has ~ 15% accuracy. 

Then it can be shown (see, e.g. Starobinsky 1988) that, in the approximation 
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outlined above, (AT/T)/ oc has the following structure: 



J^(k) = -0 o (k)^F(A;)+3zc sT G(fc)jexp 

= f(k)cos(kc s r] rec ) - g(k)sm(kc s r] rec ) 
G(k) = f(k) sm(kc s r] rec ) + g(k) cos(kc s f] rec ) 



s 



1 „ / 3Q b 



0.486 , R = -7^ = 0.413 (A15) 

3(1 + R) V^lrec 

where 0o(k) is the initial perturbation spectrum, i.e. 0o(k) = 0o(k)/c(A;), c s is the 
sound velocity in the radiation component coupled with baryons (in terms of the 
light velocity) and Rs is the Silk damping scale. Following Bond 1988, we assume 
Rs = 12 A Mpc (for the values of H Q and Q b given above). This value of Rs partly 
accounts for an increase of the Silk damping scale during recombination, and the 
final result for temperature fluctuations appears to be in a good agreement with 
numerical calculations using the exact kinetic approach (Bond & Efstathiou 1987, 
Scott & White 1994, Sugiyama 1994). The dependence of the transfer functions f(k) 
and g(k) on Qf, is weak and we neglect it. Their form in the limit Q b —* calculated 
in Sahni 1984 may be well approximated by the following analytical fits in terms of 
the variable x = kr\x = ka~ l {rio) x 309 Mpc: 

f(x) = l + 0.042x 2 x < 1.55 

= 0.8991 + 0.1302x 1.55 < x < 9.379 

= — x> 9.379, 

x+ 12.74 

31n(l + ^) 

9(?) = 7-V-1Z- (A16) 

16 ~ t ~ 400 

The oscillating terms in ( |A15[ ) are the well-known acoustic oscillations in the 
baryon-photon plasma (see, e.g., Zeldovich & Novikov 1983 and references therein). 

Substituting ( |A15|) into Eq. (|), we get after some algebra: 

x(2-K 2 B) 2 c(k 1 )c(k 2 ) e~^ +k ^ R s (F(k 1 )-3ic s ^G(k 1 )j (^F(k 2 ) - 3ic s ^G(k 2 )\ 

B 2 f vo i f°° a f°° u f 1 i 5(2A;ifc 2 + (kf + k%)u)\ „ w , , -L(t»+&u& 

= -<b L dn "l dkl l dk2 L du 2M+ iFmrkid c(kMh)e 2 
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x[F(A; 1 )F(A; 2 )ioH-3c s (F(k 1 )G(k 2 )(k 1 u + k 2 ) + G(k 1 )F(k 2 )(k 1 + k 2 u)) ^ W)in ' 



+9c 2 G'(£; 1 )G(£; 2 ) 
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' h{ w )(v - Vrec) 2 {ki + k 2 u)(k!U + k 2 ) ji(w)u\ 



w 2 w 



w = yk 2 + k 2 + 2k\k 2 u (77 — r] rec ) . (A17) 
Calculating this integral numerically, we obtain 

Iuuu « -2.2 £? 2 . (A18) 



■ T 



Hoc 



Note for comparison that if we did not take into account the increase of 
due to the Silk and Doppler effects, i.e. if we calculated the integral (|A17|) in 
the limit f(k) = l, g(k) = (but with the exact c(k)) we would get the answer 
luuu — —0.9 B 2 . So, the account of all effects increase the answer by a factor of 2.4. 



B. Sixth order contribution 



Beyond the fourth order term, the next contribution to the skewness arises 
at the sixth order in 0o- in this section, we outline the calculation for the most 
significant of the sixth order terms, ((AT/T)^\AT/T)^\AT/T)^}. Substituting 
the expression for (AT/T)[; 2) from flEjjj) and (TAD]) we obtain 



X dr]2V2 I / ^[5^(k 2 ,k 2 ) + g(k 2 ,k 2 )]e 4 

X /J dW3 / ^ / ^l5P{*>K) + Q(*i,K)W 

x (0(k 1 )0(k' 1 )0(k 2 )0(k 2 )0(k 3 )0(k' 3 )), (Bl) 

Assuming the initial potential fluctuations to be a gaussian random field (as 
predicted by most inflationary scenarios) we have 



(0(k 1 )0(k' 1 )0(k 2 )0(k' 2 )0(k 3 )0(k' 3 )) = <p 2 (ki)(j) 2 (k2)(f> 2 {k 3 ) 
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x (5(ki + k' 2 ) 5(k 2 + k 3 ) 5(k 3 + k[) + 7 permutations) (B2) 

where we ignore the self-coupling terms. We assume a scale-invariant initial 
spectrum[] with <p 2 (k) = 2ir 2 Bc 2 (k)/k 3 , where B is a normalisation constant, and 
introduce the notations u = cos#i, v = cos 9 2 and u> = cos # 3 for brevity. The 
expressions can be simplified by making the approximation 



r° VidVi f V2dv2 r r/3c/r/3e ifclfe -^ l)u e ^fc2fe -^^e ifc3li/| - i/;i ' , 

" f]rec " T)rec " Tjrec 

p7]o roo roo 

/ r/ 3 x / dx dy j^nx jk 2 { v -x)u e -ik 3 w y 

JO J— 00 J— oo 



^(2vr) 2 «y(A:iu - M <5(^ - A; 3 w) = vr 2 ^ 5(A: l2 - A; 2z ) 5(k 2z - k 3z ) (B3) 



where x = r\ 2 — rji, y — rj 3 — 771. Substituting (|B3| ) and ( p2|) into expression (pi|) , 
(73(0) can be written as (|13|). 

As in the case of the fourth order calculations, it is convenient to split C3 (0) 
into pieces involving the local and the non-local parts of & 2 \ We express 



^3^(0) — Iqqq + 3Iqqp + SIqpp + Xppp. (B4) 

Splitting the expression for cjf\o) as described by (|B"4|), we obtain the following 
results for the constituent terms: 



/ R \ 3 /> /■27T /"27T /"l /"l rl \?n-, V-^j ^.2 

ZqQQ = 4( — ) dkk 3 c 2 (k) da d/3 du dv dw c 2 (— ) c 2 (— ) — 
" " " \21/ J Jo Jo Jo Jo JO M u> w 

x cosfli cosfl 2 cos0 3 = 1.4 x lO 6 ^ 3 , (B5) 



B 



)<3 /■ /• 27T ^ 27T /* 1 /* 1 /* 1 i i 

/ dfc k 3 c 2 {k) / da / d[3 du dv dw c 2 (— ) c 2 (— ) 
J Jo Jo Jo JO JO It If 



QQP 

r 2 cos 6 1 ! — 2-uf 
r 2 — 2wu cos #i 



kv x o,kv , v 2 
uw 



cos do cos 0, = -8.1 x KT5 



5 D 3 



(B6) 



1 Our calculations can be trivially extended to non-scale-invariant spectra by incorporating the 
appropiate form for </> 2 (&). 
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fB\ 3 
1 QPP = 25 (-) 

r 2 cos 9i — 2uv 
r 2 — 2uv cos Q\ 



dk k 3 c 2 (k) 

S 2 COS #2 — 



2tt 



2tt 



2wv 



s 2 — 2wv cos 9o 



da d/3 du dv 
Jo Jo Jo 

cos#s = 1.19 x 10 7 B 2 



,kv s 



,kv s 



v 



dw c 2 (— ) c 2 ( — ) — 
u w uw 



(B7) 



Z PPP = - (™) 3 f dk k 3 c 2 (k) [ 2n da [ 27V d(3 ^ du I' dv ['dw c 2 ( — ) c 2 (— ) — 
2 V 21 / J Jo Jo Jo Jo Jo u w uw 

r 2 cos Q\ — 2uv 

X , 

r 2 — 2uv cos #i 

where we have used the following notations 

cos 6*i = uv + uvcosa, cos 9 2 = uw + uw cos /3, cos 3 = + cos /3, 
r 2 = u 2 + v 2 , s 2 = v 2 + w 2 , t 2 = w 2 + u 2 . (B9) 

We have evaluated the expressions ( p5|) , ( p6|) , ( p7|) and ( |B8|) numerically to estimate 
the contribution to the skewness at the 6th order in the potential fluctuation <p . 



s cos 62 — 2wv 
s 2 — 2wv cos 9 2 



t cos 6*3 — 2uw 
t 2 — 2uw cos 9* 



-3.68 x 10 7 B 2 



:bs) 
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